Fractional Josephson vortices at YBaL2Cu 3 7 _ x grain boundaries 
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We report numerical simulations of magnetic flux patterns in asymmetric 45° [001]-tilt grain 
boundaries in YBa2Cu307_a; superconducting films. The grain boundaries are treated as Josephson 
junctions with the critical current density j c (x) alternating along the junctions. We demonstrate 
the existence of Josephson vortices with fractional flux quanta for both periodic and random j c (x). 
A method is proposed to extract fractional vortices from experimental flux patterns. 
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Numerous recent studies of electromagnetic properties 
of grain boundaries in high-T c superconducting films are 
driven by necessity to probe the fundamental symmetry 
of the order parameter and the flux quantization [Q-^) . 
Although interpretation of the results is nontrivial, most 
of the data can be understood in terms of the conven- 
tional model of a strongly coupled Josephson junction 
P). The asymmetric 45° [001]-tilt grain boundary in 
YBa2Cu307-a; films is a notable exception of this rule. 
We show in this letter that these boundaries should have 
unusual electromagnetic properties unprecedented for the 
physics of standard Josephson junctions. 

Experimentally the asymmetric 45° [001]-tilt grain 
boundaries in YBa2Cu307_ x films exhibit an anoma- 
lous dependence of the critical current I c on the ap- 
plied magnetic field H a [§j6|- Contrary to the usual 
Fraunhoffer-type dependence with a major central peak 
at H a = and minor symmetric side-peaks, these bound- 
aries demonstrate a pattern without the central major 
peak. Instead, two symmetric major side-peaks appear 
at certain fields H a = ±H p ^0 Other remarkable 

feature is the spontaneous disordered magnetic flux gen- 
erated at the asymmetric 45° [001]-tilt grain boundaries 
in YE^CuaCv-z films Q. It is worth noting that the 
spontaneous flux is observed only in samples exhibiting 
the anomalous dependence I c {H a ). 

Clearly, the major side-peaks reveal a specific hetero- 
geneity of the Josephson properties. Indeed, a fine scale 
faceting of grain boundaries in YBCO thin films has 
been recorded by the transmission electron microscopy 
P,p|~pr| . The facets have a typical length-scale I of the or- 
der of 10-100 nm and a wide variety of orientations. This 
grain boundary structure combined with a predominant 
d x 2_ y 2 wave symmetry of the order parameter 12 lif| 
form a basis for understanding both the anomalous de- 
pendence I c (H a ) and the spontaneous flux [0j|Jl^,|l^] . 

In the case of a d x 2_ y 2 wave superconductor the phase 
difference of the order parameter across the grain bound- 
ary consists of two terms. The first term ip(x) is caused 
by a magnetic flux inside the junction and the second 
term a(x) is caused by a misalignment of the anisotropic 



banks of the junction. The Josephson current density 
j(x) depends on the total phase difference (p(x) + a(x). 
Assuming j(x) oc sin [y(x) + a(x)] one can develop a 
model of the electromagnetic properties of the grain 
boundaries in YBCO films f|. Values of the phase a(x) 
depend on the relative orientation of the neighboring 
facets. In the case of an asymmetric faceted 45° grain 
boundary we have an interchange of a = and n and 
j(x) — jc(x) sbx(p(x), where the alternating critical cur- 
rent density j c (x) oc cos a (a;). 

In this letter we report numerical simulations of flux 
patterns in the asymmetric 45° [001]-tilt grain boundaries 
in YBa2Cu307_ x superconducting films. The boundaries 
are treated as Josephson junctions with an alternating 
critical current density j c {x). We find two types of frac- 
tional Josephson vortices for each stationary state with 
a spontaneous flux in the grain boundaries, which exists 
for both periodic and random sequences of facets. One 
type of vortices contains the magnetic flux <j>\ < 4>o/2; 
the other type carries > 4>q/2 with a complementarity 
condition <f>i + (f>2 — <Pa, where </>o is the flux quantum. 
We suggest a method to extract the fractional vortices 
from the data on flux patterns. 

The alternating dependence j c {x) is imposed by a par- 
ticular sequence of facets along the boundary and there- 
fore j c {x) has the same typical length-scale / as the facets. 
It is convenient for the further analysis to write j c {x) as 



Oc) [! + <?(*)], 



(1) 



where (j c ) is the average value of the critical current den- 
sity over distances L>!: 



(jc) = — J jc{x)dx. 



(2) 



The dimensionless function g(x) in Eq. (Q) character- 
izes the Josephson properties of the grain boundary. This 
function alternates with a typical length-scale I and has 
a zero average: (g(x)) = 0. The maximum amplitude 
of alternations of \g(x)\ may Vciry from |<7(#)|max > 1 to 
|<?(a;)|max 1 depending on details of the structure of the 



1 



faceted grain boundary. We assume that A < I < Aj, 
where A is the London penetration depth and 



A 2 = 



16^ 2 A(j c > 



(3) 



is an effective Josephson penetration depth. With this 
notation, the phase difference <p(x) satisfies 



A 2 ,^" - [f + g(x)]smip = 0. 



(4) 



A model grain boundary with a periodic critical cur- 
rent density j c (x) has been considered analytically by 
means of a two-scale perturbation theory which requires 
I -c A j ]l5| ]. In this approximation, the phase (p(x) is a 
sum of a smooth part ip{x) with a length scale Aj and a 
rapidly oscillating part with a length scale I and a 
small amplitude \£(x)\ <C f: 



<p(x) = ip(x) +£(x). 

We have for the phases ip{x) an( i £( x ) : 

Aj-0" — simp + 7 sin ip cos ip = 0, 
£0*0 = Cg(^)sinV' ) 



(5) 



(6) 
(7) 



where the function £ g (x) is given by 

A$% = 9(x), (8) 
and the dimensionless parameter 7 > is defined as 

1 = -{g(x)t g (x)) = A*(Z' g 2 ). (9) 

It is worth mentioning that both £ g (x) and 7 depend only 
on the spatial distribution of j c and therefore characterize 
the individual Josephson properties of a particular grain 
boundary. We stress that this approximation is valid if 
I <C Aj and |£(a;)| -C \ip{%)\ (the latter condition results 
in \g{x)\ « 47r 2 A 2 /; 2 ). 

In the framework of the two-scale perturbation theory, 
a single Josephson vortex is described by the solution of 
Eq. (^) under the boundary conditions 0'(±oo) = 0. The 
latter can be written as smtp± (1 — 7cos0±) = 0, where 
'ijj± = ip(±oo). It is convenient for the further analysis to 
assume that ip- < ip+. 

In the case of 7 < 1 , there is only one single vortex so- 
lution, for which the phase tp(x) increases monotonically 
from ip- = to ip+ = 2tt. This solution describes the 
Josephson vortex with one flux quantum 0o . In the case 
of 7 > 1, the spatial distribution of the smooth phase 
ip describes two fractional vortices. For the first frac- 
tional vortex the phase ip(x) increases from 0_ = —ip-y 
to ip + — ip-y, where 



tp-y = arccos(I/7). 



(10) 



The difference ip+ — ip- = 2ip 1 and thus this vortex carries 
the flux (f>\ = ip~/(f)o/ir < (pa/2. For the second fractional 



vortex ip- = ip~, ip+ = 2ir — ip 7 , the phase difference 
being 2ir — 2ipj, and thus this vortex contains the flux 
4>2 = (1 — ip-y/rfcpo > </>o/2- These two fractional vortices 
are complementary: (pi + 4>2 = (pa- 
in our numerical study we solve Eq. (Q) exactly. We 
treat the stationary states as well as the relaxation to the 
stationary states using a time-dependent model Jl6| ] 



(p + aip - ip" + [1 + g(x)] sin ip = 0, 



(11) 



where a is a decay constant which we take from the in- 
terval 0.1 < a < 1. The term a<p in Eq. (|Tl| ) describes 
dissipation driving the system into one of the stable sta- 
tionary states described by the solutions of Eq. (Q). 

We begin our numerical simulations with verification 
of the results obtained by means of the two-scale ap- 
proximation for the grain boundary with periodic j c (x). 
To study the fractional vortices we start the simulations 
from a certain initial phase <pi{x) under the condition 
ipi(C)— (pi(0) = 2nn, where the boundary length C 3> Aj. 
In this case the numerical procedure converges well to a 
final stationary state. 

In Fig. [l] we show a stable stationary solution for a pair 
of fractional vortices. We compute <p(x) using the model 
g(x) = g sin(27ra/;) with g = 100 and I = O.lAj. The 
value of 7 calculated by means of Eq. <Mj is given by 



2 72 



9ll 



8^ 2 A 2 ' 



(12) 



This yields 7 » 1.27 and ip 7 « 0.66; thus <p x » 0.21 <p , 
4>2 ~ O.79 0o- As is seen in Fig. [TJ, the simulation gives 
the same value of ip 7 . The magnified insets in Fig. |l| 
demonstrate that (p(x) indeed consists of a smooth part 
superimposed with a small fast oscillating term. There- 
fore, the numerical simulations for single fractional vor- 
tices confirm the qualitative and quantitative results of 
the approximate analytic approach described above. 




20 x/Aj 

FIG. 1. The phase distribution ip(x) computed using 
g(x) = go sin(2nx/l) and 7 « 1.27. Two fractional vortices 
with 0i » 0.21 0o and 02 ~ 0.79 0o are clearly seen, the fine 
structure of <p(x) is demonstrated in the magnified insets. 
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Consider now a dilute chain of fractional vortices. Let 
a vortex with the flux <fii be situated somewhere in the 
chain. The phase ijj of this vortex changes from 27m — ?/> 7 
to 27m + "07 with an integer n. Therefore, one expects 
the phase of an adjacent vortex to start with the value 
27m + ip~f and to end up with 27r(n + 1) — ip 7 , the total 
phase accumulation of these two vortices being 27r. In 
other words, the chain consists of a sequence of pairs of 
vortices with fluxes tf>\ and <f>2- This qualitative picture 
is confirmed by numerically solving Eq. ([!]) . Figure |^(a) 
shows the result of such a calculation for which we took 
g(x) = 150 sin(20 7ra;/Aj), that corresponds to 7 « 2.85, 
■0 7 ~ 1.21, and the fluxes <f>i » O.390o, 4>2 ~ 0.61</>o- The 
final stationary state of our numerical procedure simulat- 
ing the relaxation process, depends on the choice of the 
initial phase tpi(x). By taking a proper non-monotonic 
dependence ipi(x), we may end up with a stationary so- 
lution shown in Fig. ||(b). A remarkable feature of this 
solution is the existence of fractional vortex-antivortex 
pairs clearly seen in the simulation of Fig. |2|(b), the pair 
with the fluxes ±0i is followed by the pair with 
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FIG. 2. Two chains of fractional vortices in a grain bound- 
ary with a periodically alternating critical current density: 
(a) an "ideal" chain, (6) a chain with vortex-antivortex "de- 
fects". Empty triangles mark the positions of the fractional 
vortices with the fluxes <j>i < </>o/2, full triangles correspond 
to cj>2 > 0o/2. The up-down orientation of triangles indicate 
the field direction of vortices. For this particular calculation 
we use g(x) = go sin(2ivx/l), go — 150, I — 0.1 Aj, which 
result in 7 ~ 2.85, i/> 7 ~ 1.21, and 4>\ ~ 0.39</>o, <j>2 ~ O.610o- 

Next, we study flux patterns for a more realistic case of 
a grain boundary with 2A^ facets and a non-periodic al- 
ternating critical current density j c {x). We treat this case 
numerically and use a stepwise g(x) defined as: g{x) = go 
if 0^ < x < bi, and g(x) = —go if bi < x < a^+i 
(i = l,...,N). It is convenient to introduce the ran- 
dom distances d~i and bi with {di) = (bi) = and a stan- 
dard deviation o~i such that bi — ai — 0.5(Z + di) and 
flj+i — bi = 0.5(7 + bi). The simulations start with an ini- 
tial phase <fi(x) matching the condition (fi(C) — (fi(0) = 
2ixn. Then the numerical procedure of solving Eq. ( |Tl| ) 
converges well to a stationary state which depends on 



both ifi (x) and a due to the flux pinning induced by the 
non-uniformity of the critical current density. 

A special role in the description of Josephson bound- 
aries with random alternating j c {x) belongs to the sta- 
tionary state if s (x) which corresponds to the zero total 
spontaneous flux and to the absolute minimum of the 
Josephson energy £{(p(x)} (defined in a standard way 
|i"7f ). Our simulations show that Lp s {x) is unique for a 
given boundary, stable, and independent either of initial 
guesses ipi(x) or of the damping constant a. Therefore, 
the phase <p s {x) can serve as a signature of each indi- 
vidual boundary. It is convenient to represent ip s (x) as 
^p s {x) = ipj + £, s ( x ) with ip 1 = const and the variable 
part £, s ( x ) having zero average, (^ s (x)) — 0, and a typi- 
cal amplitude |£ s (aOI < 7r /2- 

An example of a computed ip s (x) is shown in Fig. [I](a). 
For this simulation we took ifi (x) = const + £j (x) with 
an arbitrary small As stated above, the resulting 

phase tp s (x) is independent of (fi{x). It is worth men- 
tioning that the spontaneous self-generated flux 4> s {x) = 
4>q£,s{x) /27r has a wide range of length-scales imposed by 
the random j c {x). Note also that randomness of j c (x) re- 
sults in a considerably higher amplitude of the flux vari- 
ation 4> s (x) as compared to a periodic j c ( x )] this is seen 
from comparison of Fig. ||(a) with the insets in Fig. |[ 




(a) 
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FIG. 3. Two stationary solutions <fi(x) developed in a zero 
magnetic field for a stepwise randomly alternating g(x) for 
two initial conditions: (a) preventing, (b) stimulating creation 
of vortices (go = 200, I = 0.1 Aj, oi » 0.06/, ^ 7 « 1.515, 
<t>i « O.480o, (t>2 « O.520 O ). 

It follows from Eq. that the stationary solution 
ip s {x) generates two series of solutions having the same 
Josephson energy as <p s (x): <fin{x) — 2im + <p s (x) and 
fni x ) = — ip s (x), where n is an integer. The aver- 
age values (99+ (x)) = 2im + 7 and (ip~ (x)) = 2irn — ip 7 
interchange being separated by (fn( x )) ~ (Vn ( x )) = ^7 
or by (<Pn+i( x )) ~ (Pni 30 )) = 2ir — 2ip 7 , as is shown in 
Fig. U for a periodic j c { x )- The gaps between the average 
values of the stationary phases fn( x ) and Vn{ x ) allow 
for fractional vortices with the fluxes <pi = fai/j^/ir and 
4>2 = <Po ~ 4>i as solutions of Eq. (^) varying with a typ- 
ical length-scale of A,/. An example of a computed ip(x) 
with two clearly pronounced fractional vortex-antivortex 
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pairs and small varying part £(x) is shown in Fig. ||(b). 

Consider now a typical flux (phase) pattern for a chain 
of vortices for a randomly alternating j c (x) . Assume that 
the chain starts with a domain with ip(x) = <p s (x), i.e., 
the phase varies slightly (|£ s (a;)| < tt/2) around its av- 
erage value tj)j. Therefore, ipy is the value of (tp(x)) at 
the "tail" of the neighboring vortex or antivortex. If the 
neighbor carries the flux 02, the average phase should 
increase from -0 7 to 2?t — in the neighbor's domain. 
Alternatively the neighbor may carry the flux —<j>\ gen- 
erating a decrease of (ip(x)) from ip^ to —ipy In general, 
this flux pattern is similar to the one arising for a pe- 
riodic j c (x). However, the relatively large amplitudes 
of the spontaneous flux and the absence of periodicity 
within the chain may mask the fractional vortices. 
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one of these lines everywhere, except a few relatively 
sharp jumps from one line to the next; in particular, the 
"signature" ip s (x) is nested at ip-y'- <Ps(x) — ^P-/ + £,s{x). 
The jumps (or vortices) should be centered at tp(x) = mi 
[where according to Eq. ([!]) ip"(x) = 0]. Then we take 
a domain situated between lines im and ir(n + 1) and 
form ip v (x) — ip(x) =F £ s (%)> choosing the minuses if the 
random parts of <p(x) and ip s (x) are identical, and the 
plus otherwise. The curve ip v (x) shown by a thick curve 
at Fig. ^(b) is remarkably smooth; clearly it represents 
the fractional vortices described above within the model 
with periodic j c {x). 

In conclusion, we have shown by numerical simula- 
tions that two types of fractional vortices with fluxes 
0i < o /2 and 2 = 0o ~ <t>i > 71 7 2 exist at 45° [001]-tilt 
grain boundaries in YBa2Cu307_a; films exhibiting spon- 
taneous flux in zero field cooled samples. Faceted grain 
boundaries are treated as Josephson junctions with alter- 
nating critical current density. We show how to extract 
fractional vortices from the data on spontaneous flux dis- 
tribution. 

One of us (RGM) is grateful to J.R. Clem, V.G. Kogan, 
and J. Mannhart for useful and stimulating discussions. 
This research is supported in part by grant No. 96-00048 
from the United States - Israel Binational Science Foun- 
dation (BSF), Jerusalem, Israel. 



x/Aj 



FIG. 4. (a) The "signature" y> B (x) and the phase f(x) 
with the total flux = 2 0o at the boundary calculated for 
go = 90, / = 0.1 Aj, oi ~ 0.015 Z, which gives ip-, m 1.21, and 
0i ~ O.390o, 02 ~ 0.6100- (6) The thin line depicts <p(x), the 
thick line depicts the phase tp v (x) generated by the fractional 
vortices and extracted from the phase ip(x). 

At the top panel of Fig. || we show (p(x) obtained nu- 
merically for go = 90, I = 0.1 A,/, and 07 = 0.015 1, the 
calculated value of -0 7 is 1.21. The bottom curve is the 
"signature" (p s (x) corresponding to the zero total flux, 
= 0. The upper curve is calculated with the same set 
of input parameters for = 20o. Even a brief examina- 
tion of the curves shows a striking correlation between 
the two: there are domains (e.g., —10 < x/Aj < 10) in 
which the "noise" patterns are nearly identical, whereas 
in others (e.g., 10 < x/Aj < 30) the patterns repeat each 
other being flipped. This suggests that one can extract 
the smooth part of the upper curve f(x) by properly 
subtracting the "signature" ip s (x). 

The subtraction is done as follows. First, we draw the 
straight lines 27rn±?/' 7 at the graph of (f(x), see Fig. ||(b). 
We see that "random" variations of f(x) are nested on 
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